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We  report  the  discovery  in  Lake  Cuitzeo  in  central  Mexico  of  a 
black,  carbon-rich,  lacustrine  layer,  containing  nanodiamonds,  mi¬ 
crospherules,  and  other  unusual  materials  that  date  to  the  early 
Younger  Dryas  and  are  interpreted  to  result  from  an  extraterres¬ 
trial  impact.  These  proxies  were  found  in  a  27-m-long  core  as  part 
of  an  interdisciplinary  effort  to  extract  a  paleoclimate  record  back 
through  the  previous  interglacial.  Our  attention  focused  early  on 
an  anomalous,  10-cm-thick,  carbon-rich  layer  at  a  depth  of  2.8  m 
that  dates  to  12.9  ka  and  coincides  with  a  suite  of  anomalous 
coeval  environmental  and  biotic  changes  independently  recog¬ 
nized  in  other  regional  lake  sequences.  Collectively,  these  changes 
have  produced  the  most  distinctive  boundary  layer  in  the  late 
Quaternary  record.  This  layer  contains  a  diverse,  abundant  assem¬ 
blage  of  impact-related  markers,  including  nanodiamonds,  carbon 
spherules,  and  magnetic  spherules  with  rapid  melting/quenching 
textures,  all  reaching  synchronous  peaks  immediately  beneath  a 
layer  containing  the  largest  peak  of  charcoal  in  the  core.  Analyses 
by  multiple  methods  demonstrate  the  presence  of  three  allotropes 
of  nanodiamond:  n-diamond,  i-carbon,  and  hexagonal  nanodia¬ 
mond  (lonsdaleite),  in  order  of  estimated  relative  abundance.  This 
nanodiamond-rich  layer  is  consistent  with  the  Younger  Dryas 
boundary  layer  found  at  numerous  sites  across  North  America, 
Greenland,  and  Western  Europe.  We  have  examined  multiple  hy¬ 
potheses  to  account  for  these  observations  and  find  the  evidence 
cannot  be  explained  by  any  known  terrestrial  mechanism.  It  is, 
however,  consistent  with  the  Younger  Dryas  boundary  impact 
hypothesis  postulating  a  major  extraterrestrial  impact  involving 
multiple  airburst(s)  and  and/or  ground  impact(s)  at  12.9  ka. 

black  mat  |  cosmic  impact 

We  present  data  from  Lake  Cuitzeo  in  central  Mexico  (19.94  °N, 
101.14  °W)  in  support  of  evidence  for  the  Younger  Dryas 
(YD)  impact  hypothesis,  as  first  presented  at  the  2007  Meeting 
of  the  American  Geophysical  Union  in  Acapulco,  Mexico. 
There,  a  consortium  of  scientists  reported  geochemical  and 
mineralogical  evidence  from  multiple  terrestrial  sites  ascribed 
to  extraterrestrial  (ET)  impacts  and/or  airbursts  (1).  Their  first 
evidence  was  the  discovery  at  well-dated  Clovis-era  archeological 
sites  in  North  America  of  abundant  magnetic  spherules  (MSp) 
and  carbon  spherules  (CSp)  in  a  thin  layer  (0.5  to  5  cm)  called 
the  Younger  Dryas  boundary  layer  (YDB),  dating  to  12.9  ± 
0.1  kaBP  (calibrated,  or  calendar  years)  or  10.9  14CkaBP 
(radiocarbon  years)1  (1-3).  The  YDB  is  commonly  located  di¬ 
rectly  beneath  or  at  the  base  of  an  organic-rich  layer,  or  “black 
mat,”  broadly  distributed  across  North  America  (1).  Later,  abun¬ 
dant  nanodiamonds  (NDs)  were  discovered  by  Kennett  et  al. 
(2,  3)  in  the  YDB  layer  at  numerous  locations.  NDs  also  were 


detected  at  the  margin  of  the  Greenland  Ice  Sheet  in  a  layer 
that  dates  to  the  approximate  YD  onset  (4).  These  discoveries 
led  to  the  hypothesis  that  one  or  more  fragments  of  a  comet  or 
asteroid  impacted  the  Laurentide  Ice  Sheet  and/or  created  at¬ 
mospheric  airbursts  (1)  that  initiated  the  abrupt  YD  cooling 
at  12.9  ka,  caused  widespread  biomass  burning,  and  contributed 
to  the  extinction  of  Late  Pleistocene  megafauna  and  to  major 
declines  in  human  populations  (5). 

Some  independent  workers  have  been  unable  to  reproduce 
earlier  YDB  results  for  MSp,  CSp,  and  NDs  (6-8),  as  summarized 
in  a  “News  Focus”  piece  in  Science  (9),  which  claims  that  the 
YDB  evidence  is  “not  reproducible”  by  independent  researchers. 
Refuting  this  view,  multiple  groups  have  confirmed  the  presence 
of  abundant  YDB  markers,  although  sometimes  proposing  alter¬ 
nate  hypotheses  for  their  origin.  For  example,  Mahaney  et  al. 
(10-12)  independently  identified  glassy  spherules,  CSps,  high- 
temperature  melt-rocks,  shocked  quartz,  and  a  YDB  black  mat 
analogue  in  the  Venezuelan  Andes.  Those  authors  conclude 
the  cause  was  “either  an  asteroid  or  comet  event  that  reached 
far  into  South  America”  at  12.9  ka.  At  Murray  Springs,  Arizona, 
Haynes  et  al.  (13)  observed  highly  elevated  concentrations  of 
YDB  MSp  and  iridium.  Abundances  of  MSp  were  340  x  higher 
than  reported  by  Firestone  et  al.  (1)  and  iridium  was  34  x  higher, 
an  extraordinary  enrichment  of  3,000  x  crustal  abundance.  Those 
authors  stated  that  their  findings  are  “consistent  with  their  (Fire¬ 
stone  et  al.’s)  data.”  In  YDB  sediments  from  North  America  and 
Europe,  Andronikov  et  al.  (2011)  reported  anomalous  enrich¬ 
ments  in  rare  earth  elements  (REE)  and  “overall  higher  concen¬ 
trations  of  both  Os  and  Ir  [osmium  and  iridium]”  that  could 
“support  the  hypothesis  that  an  impact  occurred  shortly  before 
the  beginning  of  the  YD  cooling  12.9  ka.’4.  Tian  et  al.  (14) 
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observed  abundant  cubic  NDs  at  Lommel,  Belgium,  and  con¬ 
cluded  that  “our  findings  confirm  ...  the  existence  of  diamond 
nanoparticles  also  in  this  European  YDB  layer.”  The  NDs  occur 
within  the  same  layer  in  which  Firestone  et  al.  (1)  found  impact- 
related  materials.  Similarly,  at  a  YDB  site  in  the  Netherlands, 
Van  Hoesel  et  al.§observed  “carbon  aggregates  [consistent  with] 
nanodiamond.”  Recently,  Higgins  et  aV  independently  an¬ 
nounced  a  4-  to  4.5-km-wide  YDB  candidate  crater  named  Cor- 
ossol  in  the  Gulf  of  St.  Lawrence,  containing  basal  sedimentary 
fill  dating  to  12.9  ka.  If  confirmed,  it  will  be  the  largest  known 
crater  in  North  and  South  America  within  the  last  35  million  years 
Because  of  the  controversial  nature  of  the  YD  impact  debate, 
we  have  examined  a  diverse  assemblage  of  YDB  markers  at  Lake 
Cuitzeo  using  a  more  comprehensive  array  of  analytical  techniques 
than  in  previous  investigations.  In  addition,  different  researchers  at 
multiple  institutions  confirmed  the  key  results. 

Lake  Cuitzeo 

Covering  300-400  km2,  Cuitzeo  is  the  second  largest  lake  in 
Mexico,  located  at  high  elevation  of  1,820  m  in  the  northern  part 
of  the  state  of  Michoacan  within  the  tectonically  active  Trans 
Mexico  volcanic  arc  (15)  (SI  Appendix,  Fig.  1).  Situated  in  the 
tropics,  this  lake  currently  experiences  a  semiarid  climate  with 
annual  temperatures  ranging  from  10  to  28  °C  (avg  19  °C)  and 
annual  rainfall  ranging  from  60  to  100  cm.  This  shallow  lake  cur¬ 
rently  varies  in  depth  from  0.8  to  2.2  m  (avg  1.9  m). 

Results 

Sedimentary  Sequence.  A  27-m-long,  10-cm-diameter  core  was  ex¬ 
tracted  in  1997  from  thick  deposits  in  Lake  Cuitzeo  as  part  of  an 
interdisciplinary,  multiproxy  effort  to  acquire  a  detailed  paleocli- 
mate  record  extending  back  to  the  last  interglacial  [130,000 
(130  kyr)]  (15).  The  core  consists  of  interbedded  sands,  silts,  clays, 
and  epiclastites,  along  with  6-kyr-old,  20-cm-thick  tephra  between 
1.7  and  1.4  m  and  31.5-kyr-old,  20-cm-thick  tephra  between  4.7  and 
4.5  m,  with  several  more  volcanic  deposits  below  10  m.  A  conspic¬ 
uous,  dark,  carbon-rich  layer,  dominantly  comprised  of  clay  and 
silt,  occurs  between  2.82  and  2.50  m  (Fig.  1)  and  is  the  focus  of 
this  study  because  of  its  similarity  to  the  black  mat  at  YDB  sites 
across  North  America.  Sediment  samples  of  approximately  1  cm 
thickness  were  taken  every  5  cm  across  the  critical  section  between 
2.80  and  2.65  m  and  at  10  cm  intervals  above  and  below  this  sec¬ 
tion.  These  samples  were  quantitatively  analyzed  for  diatoms  and 
pollen  assemblages,  carbonate  (%TIC),  organic  carbon  (%TOC), 
bulk  major-element  composition,  stable  carbon  isotopes  (both  or¬ 
ganic  and  inorganic),  organic  nitrogen,  MSp,  NDs,  CSp,  charcoal, 
and  aciniform  soot. 

Chronology.  Previously,  Israde  et  al.  (15)  published  an  age-depth 
model  for  the  uppermost  9  m  at  Lake  Cuitzeo  comprised  of  16 
accelerator  mass  spectroscopy  (AMS)  14  C  dates  on  bulk  sediment 
and  used  in  a  linear  interpolation  with  the  YD  onset  identified  at 
approximately  2.8  m.  To  test  and  refine  that  model,  we  acquired 
six  more  AMS  14  C  dates  on  bulk  sediment,  for  a  total  of  22  dates, 
and  calibrated  them  using  the  IntCal04  calibration  curves  in  Cal- 
Pal0711  (Fig.  1;  SI  Appendix,  Table  1).  A  20-cm-thick  tephra  layer 
at  4.7  to  4.5  m  has  been  identified  as  the  Cieneguillas  rhyolitic 
tephra,  originating  from  nearby  Las  Azufres  volcano  and  14  C  da¬ 
ted  by  others  at  three  locations  to  approximately  31  ka 
(26.8  ±  0.9  14C  ka)  (16).  The  age  of  this  tephra  serves  as  an 
anchor  for  the  chronology.  The  dates  from  9  to  3.35  m  and  from 
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Fig.  1.  (Left).  Lake  Cuitzeo  lithostratigraphy  from  4.0  to  2.0  m.  Red  brackets 
indicate  the  carbon-rich  layer  corresponding  to  the  YD.  Blue  tick  marks  at  left 
indicate  sample  depths.  (Right)  Graph  of  calibrated  14C  dates.  A  regression 
polynomial  (black  line)  of  accepted  dates  (red  circles)  and  tephra  date  (black 
dot);  blue  circles  are  excluded  dates.  Error  bars  are  less  than  circle  widths. 
Dark  gray  band  denotes  YD  interval;  lighter  gray  band  corresponds  to  inter¬ 
val  between  4.0  and  2.0  m.  Cal  ka  BP,  calibrated  kiloannum  before  present; 
char,  charcoal. 

2  to  0  m  show  a  relatively  consistent  linear  increase  with  depth, 
ranging  in  age  from  approximately  46  to  0  kyr. 

The  samples  at  3.35  and  1.95  m  have  calibrated  ages  of  18.8 
and  9.9  ka,  respectively,  consistent  with  the  linear  extrapolation 
of  the  rest  of  the  core.  However,  six  samples  between  these  two 
levels  provided  radiocarbon  ages  older  than  the  interpolation 
predicts.  They  represent  a  major  radiocarbon  reversal  of  thou¬ 
sands  of  years,  with  older  sediment  overlying  younger,  a  situation 
that  can  result  from  reworking  of  older  organic  material.  The  re¬ 
versal  begins  with  a  date  of  18.8  ka,  shifts  anomalously  older  by 
approximately  20  kyr  above,  and  then  normalizes  to  9.9  ka  higher 
in  the  section.  At  2.75  m  within  this  interval,  total  organic  carbon 
(TOC)  is  15.8  wt%,  the  highest  percentage  in  approximately 
100  kyr  (15).  This  sample  yields  a  date  of  approximately  32  ka, 
but  linear  interpolation  indicates  it  should  date  to  approximately 
13  ka,  a  difference  of  approximately  20  kyr.  Accounting  for  this 
shift  requires  major  contamination  of  the  TOC  by  radiocarbon- 
dead  or  very  old  carbon  (92  wt%).  Currently,  the  source  of  this 
old  carbon  remains  unclear. 

To  compensate  for  these  anomalously  old  radiocarbon  dates, 
we  excluded  the  six  dates  that  form  the  reversal  between  3.10  and 
2.05  m  and  utilized  the  remaining  16  dates  to  generate  an  age- 
depth  curve  with  a  fifth-order  polynomial  regression*  .  The  result¬ 
ing  curve  (Fig.  1),  predicts  that  the  12.9-ka  YD  onset  is  at  a  depth 
of  approximately  2.9  to  2.7  m,  consistent  with  the  earlier  identi¬ 
fication  by  Israde  et  al.  (15). 

Biostratigraphy.  Although  the  stratigraphic  position  of  the  YD  on¬ 
set  has  been  reasonably  extrapolated  using  numerous  14  C  dates, 
as  is  standard  practice,  we  have  examined  other  stratigraphic  data 
to  assist  with  this  placement  of  the  YD  onset.  This  is  achieved 
using  biostratigraphic  correlation  of  pollen  sequences  from  Lake 
Cuitzeo  with  those  from  YD -aged  regional  lakes  that  have  been 
independently  dated.  Islebe  and  Hooghiemstra  (17)  reported 
that  evidence  for  YD  climate  change  is  either  present  or  likely 
in  some,  but  not  all,  lakes  in  Mexico,  Guatemala,  Costa  Rica, 
Colombia,  Ecuador,  and  Peru.  Of  those  mentioned,  we  have 
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examined  pollen  records  from  Lake  Peten  Itza  in  Guatemala  (18, 
19),  La  Chonta  Bog  in  Costa  Rica  (20),  and  Lake  La  Yeguada  in 
Panama  (21)  (SI  Appendix ,  Fig.  2).  The  pollen  sequences  from 
these  lakes  reflect  climate  changes,  including  the  Late  Glacial, 
Bplling-Allerpd  (BA),  (14.5  to  12.9  ka),  YD  stadial  (12.9  to 
11.5  ka),  and  the  Holocene  (11.5  ka  to  present),  comprising  a 
distinctive  cold-warm-cooler-warmer  climatic  sequence.  The  YD 
interval  was  specifically  identified  by  previous  investigators  at 
Lakes  La  Chonta  and  Peten  Itza. 

Pollen  changes  in  Lake  Peten  Itza  at  the  YD  onset  display  a 
decrease  in  Quercus  (oak),  a  persistence  of  Pinus  (pine),  a  dom¬ 
inance  of  Poaceae  (grasses),  and  low  diversity  and  productivity 
for  plants  in  general  (19).  At  or  near  the  YD  onset,  the  record 
at  Peten  Itza  exhibits  exceedingly  large,  abrupt,  and  unprece¬ 
dented  changes  (both  in  magnitude  and  rate  of  change)  for  tem¬ 
perature,  rainfall,  and  biotic  turnover.  These  changes  produced 
the  most  distinctive  layer  in  the  Late  Quaternary  record  (18,  19). 

For  Lake  La  Yeguada,  although  the  YD  episode  was  not 
expressly  identified,  workers  there  recognized  a  major,  abrupt 
environmental  and  ecological  change  (a  “time  of  crisis”)  close  to 
the  onset  of  the  YD  at  approximately  12.8  ka  (10.8  14C  ka)  (20). 
This  is  reflected  in  dramatic  changes  in  the  pollen  and  diatom 
records,  biotic  turnover,  clay  mineralogy,  sedimentary  geochem¬ 
istry,  and  particulate  carbon  flux.  At  the  same  time,  Quercus  and 
Myrtaceae  (myrtle)  were  replaced  by  Poaceae.  These  changes  re¬ 
present  the  most  distinctive  layer  in  this  record  (21). 

In  these  regional  lakes,  the  pollen  records  typically  form  a 
“peak-trough-peak”  pollen  pattern  with  the  trough  representing 
very  low  pollen  levels  during  the  YD  (Fig.  2;  SI  Appendix , 
Table  2).  At  Cuitzeo,  total  pollen  and  Quercus  reveal  a  similar 
pattern  in  which  high  pollen  abundances  below  2.9  m  correspond 
to  the  BA  when  warmer  temperatures  supported  abundant  biota 
around  the  lake.  Low  pollen  abundances  in  the  trough  corre¬ 
spond  to  the  YD  interval.  For  all  lakes,  the  unique,  distinctive 
BA  pollen  peaks  are  among  the  largest  in  the  last  40  to  100  kyr 
and  are  followed  during  the  YD  by  some  of  the  lowest  total  pollen 
values,  consistent  with  a  cooler  and/or  drier  climate.  After  the 
YD,  a  rebound  of  varying  magnitude  occurred  at  all  lakes  after 
the  Holocene  began  at  11.5  ka. 

Because  climate  change  also  affects  diatom  populations,  we 
compared  the  diatom  record  at  Cuitzeo  with  other  lakes,  and  as 
with  pollen,  the  diatom  record  reveals  a  period  of  extraordinary 
change  at  the  YD  onset.  Stephanodiscus  niagarae  andAulacoseira 
spp.  display  major  YD  abundance  peaks  that  are  among  the  lar¬ 
gest  in  the  last  100  kyr  (SI  Appendix,  Fig.  3  and  Biostratigraphy). 
For  Cuitzeo  diatom  assemblages,  we  also  plotted  the  change  in 
diversity  (“5  diversity”),  defined  as  the  value  derived  from  sub¬ 
tracting  the  total  number  of  diatom  species  in  one  sample  from 
the  adjacent  sample  above  it.  This  demonstrates  that  the  greatest 
change  in  diversity  within  >50  kyr  occurred  at  the  YD  onset  at 
2.8  m  in  Lake  Cuitzeo  (5  diversity  in  SI  Appendix,  Fig.  3),  and  this 
correlates  well  with  the  Lake  La  Yeguada  record  where  the  great¬ 
est  turnover  in  diatom  species  occurred  at  approximately  12.8  ka 
(10.8  14C  ka)  (21).  In  addition,  for  Lakes  Cuitzeo  and  La  Yegua¬ 


da,  especially  large  YD  peaks  are  evident  in  other  aquatic  taxa, 
including  cattails  (Typha)  and  the  algal  forms  Botryococcus  and 
Coelastrum  (SI Appendix,  Fig.  4).  High  abundances  for  these  taxa 
are  consistent  with  major  ecological  change  at  the  YD  onset. 

In  Fig.  2,  we  compare  the  pollen  records  to  a  temperature 
proxy  (5lsO)  from  a  Greenland  ice  core,  GISP2  (22).  We  have 
also  correlated  paleoceanographic  records  from  the  Cariaco 
Basin  in  the  Caribbean  (Fig.  2),  in  which  titanium  represents  ter¬ 
rigenous  input  due  to  continental  runoff  (23)  and  molybdenum 
varies  in  response  to  deglacial  climate  change  (24).  The  YD  onset 
is  identified  in  all  three  records  and  corresponds  well  with  the 
pollen  records  from  Cuitzeo  and  other  lakes. 

In  summary,  from  widely  separated  lakes  in  the  highlands  of 
Costa  Rica  to  the  lowlands  of  Guatemala  and  Panama,  there 
is  only  one  stratigraphic  interval  that  displays  extraordinary  en¬ 
vironmental  and  biotic  changes,  and  in  each  case,  this  interval 
occurs  at  or  near  the  YD  onset.  For  Lake  Cuitzeo,  the  age-depth 
model  indicates  the  YD  onset  occurs  between  1.95  and  3.35  m, 
representing  a  9-kyr  span.  Within  this  span,  only  one  level  dis¬ 
plays  extraordinary  environmental  and  biotic  changes,  as  in  other 
regional  lakes,  and  that  level  is  at  2.8  m.  Therefore,  we  conclude 
that  the  Cuitzeo  age-depth  model  is  robust  and  that  the  YD  onset 
is  correctly  identified  at  2.8  m. 

Sedimentary  Geochemistry.  Our  attention  was  first  drawn  to  the 
anomalous  interval  of  unusually  high  values  of  TOC  (5-16%)  oc¬ 
curring  between  4.0  and  2.6  m,  particularly  a  TOC  of  15.8%  in  the 
thin,  1-cm-thick  layer  at  2.75  m,  close  to  the  YD  onset  (Fig.  3). 
This  value  is  the  highest  TOC  in  the  entire  27-m  core,  which  has  a 
background  average  of  only  1.2%.  We  performed  Sl3C  analyses, 
and  below  2.75  m,  there  are  minor  fluctuations  near  an  average  of 
— 2%e,  which  is  typical  for  algal  matter.  Above  2.75  m,  8UC  va¬ 
lues  increase  10  x  to  —19 %c  at  2.7  m  in  the  dark  layer,  followed 
by  heavier  values  above  2.0  m  (after  approximately  10  ka)  in  the 
Holocene.  The  2.75-m  layer  is  depleted  in  phosphorus,  producing 
a  distinct  carbon/phosphorus  ratio  (C/P)  peak  that  is  the  highest 
in  core.  This  material  may  be  analogous  to  the  YD  carbon-rich 
black  mat  observed  at  many  North  American  sites  (25).  When 
viewed  with  the  scanning  electron  microscope  (SEM),  the 
2.75-m  layer  contains  thin  millimeter-sized  interbeds  of  black  or¬ 
ganic  carbon  that  appear  without  form  or  structure.  Analysis  by 
energy-dispersive  X-ray  spectroscopy  (EDS)  indicates  these 
bands  are  almost  pure  elemental  carbon. 

Organic  matter  in  the  anomalous  interval  is  enigmatic  and  not 
the  normal  plant-derived  kerogenous  organic  matter,  dominating 
the  rest  of  the  27-m  core  over  the  last  100  kyr;  instead,  it  appears 
to  be  very  old  and  radiocarbon-dead.  Pyrolysis  analysis  (Rock/ 
Eval)  was  conducted  on  samples  from  the  carbon-rich  layer  be¬ 
tween  2.90  and  2.55  m  and  then  compared  with  a  carbon-rich 
sample  from  the  YD -aged  black  mat  at  Murray  Springs,  Arizona 
(SI  Appendix,  Table  3),  as  previously  analyzed  by  Bischoff  in 
Haynes  (25).  This  comparison  suggests  that  much  of  the  TOC 
is  unreactive  carbon,  whereas,  according  to  GC-MS  analysis, 
the  remaining  extractable  carbon  fraction  is  typical  of  immature 
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Fig.  2.  Graphs  for  pollen,  Cariaco  Basin  proxies,  and  GISP2 
temperatures.  ( A-Q  compare  Lake  Cuitzeo  pollen  abun¬ 
dances  to  two  regional  lakes.  Warmer  Bolling-Allerod 
(BA)  in  light  gray,  and  YD  in  dark  gray.  D  and  E  show  that 
Lake  Cuitzeo  Quercus  abundances  are  similar  to  those  of 
Lake  Peten  Itza.  All  lake  plots  correspond  well  to  graph  F 
from  a  Cariaco  Basin  core  displaying  ppm  abundances  of 
titanium  (orange;  smoothed  30  x)  and  molybdenum  (red; 
smoothed  3  x)  (23,  24).  Graph  G  is  a  GISP2  temperature 
proxy  plot  (%o  <5180;  smoothed  10  x  )  (22).  Black  diamonds 
are  depth  of  14C  dates.  All  graphs  are  similar,  demonstrating 
that  the  YD  onset  is  consistent  at  all  sites. 
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Fig.  3.  Changes  in  carbon  for  the  upper  6  m  of  the  Lake  Cuitzeo  sequence. 
There  YD  onset  peaks  in  TOC  wt%,  C/P,  and  <513C.  The  dark  gray  band  denotes 
the  YD  interval,  and  the  light  gray  band  is  the  interval  between  4.0  and  2.0  m. 


plant-derived  compounds,  mostly  n-alkanes.  Contamination  of 
carbon  from  petroleum  seeps,  such  as  those  in  Lake  Chapala 
300  km  to  the  west,  was  explored  as  a  possible  source  of  the 
anomalous  carbon  (26).  However,  pyrolysis  analysis  showed  no 
detectable  petroleum  hydrocarbons  in  the  sequence  analyzed. 
This  carbon  does  not  appear  to  be  derived  from  typical  immature 
plant  compounds,  and  its  origin  is  unknown. 


Fig.  5.  SEM  images  of  magnetic  impact  spherules.  ( A-B )  Magnetic  impact 
spherules  with  dendritic  surface  pattern.  (C)  Framboidal  pyrite  spherule. 
(D)  Collisional  magnetic  impact  spherules.  (£)  Light  micrograph  of  same  mag¬ 
netic  impact  spherules.  (F)  Teardrop-shaped  spherule  with  dendritic  pattern. 
(G)  Photomicrograph  of  same  MSps.  For  labels  such  as  "2.80  #3,"  "2.80"  re¬ 
presents  depth  of  sample  in  meters  and  "#3"  is  the  magnetic  impact  spherule 
number  as  listed  in  SI  Appendix,  Table  4. 


Impact  Proxies 

CSpf  Charcoal,  and  Aciniform  Soot.  Black  CSp,  20  to  260  pm,  aver¬ 
aging  90  pm  diameter,  were  observed  at  Cuitzeo,  appearing  as 
ovoid-to-round  with  cracked,  roughened  surfaces  and  typically  re¬ 
vealing  a  thin  rind,  with  spongy,  vesicular  interiors  surrounded  by 
a  smooth,  homogeneous  matrix  (SI  Appendix,  Fig.  5).  SEM-EDS 
indicates  that  CSp  are  dominantly  carbon  (>87%)  with  minor 
particulates,  such  as  Si,  Al,  and  Fe,  concentrated  in  the  rind.  They 
reach  a  significant  peak  of  approximately  680/kg  at  2.75  m  within 
the  YDB  layer.  Above  2.75  m,  the  CSp  persist  at  an  average  of 
approximately  120/kg.  No  CSp  were  detected  below  2.8  m. 

Charcoal  microparticles  (>125  pm)  were  counted  between 
3.6  and  2.2  m,  an  interval  dating  from  approximately  21  to 
10  ka  (Fig.  4),  displaying  background  levels  of  approximately 
5,000  particles/kg.  There  is  a  minor  charcoal  peak  at  3.1  m  of  ap¬ 
proximately  29,000  particles/kg,  dating  to  approximately  16  ka, 
and  there  is  a  major  rise  in  charcoal  beginning  near  the  YD  onset 
and  reaching  a  maximum  peak  at  2.65  m  of  77,000  particles/kg 
(15  x  background)  in  a  major  episode  in  biomass  burning.  The 
main  charcoal  peak  is  about  5  cm  above  the  impact  proxies  dis¬ 
cussed  below,  and  the  lack  of  tephra  within  this  interval  indicates 
the  biomass  burning  is  unrelated  to  volcanism.  There  is  also  a  ma¬ 
jor  peak  in  particulate  carbon  (charcoal)  in  Lake  Le  Yeguada  that 
dates  close  to  12.8  ka  (21),  near  the  YD  onset. 
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Fig.  4.  Markers  over  the  interval  between  3.6  and  2.2  m.  The  YD  episode 
(12.9  to  11.5  ka)  is  represented  by  dark  band.  YDB  layer  is  at  2.8  m.  NDs 
and  magnetic  impact  spherules  both  peak  at  the  YD  onset,  whereas  framboi¬ 
dal  spherules,  CSps,  and  charcoal  peak  higher  in  the  sequence.  Magnetic 
grains  peak  just  prior  to  the  YD  onset.  NDs  are  in  ppb;  Msps,  framboidal 
spherules,  CSps,  and  charcoal  are  in  no./kg;  magnetic  grains  in  g/kg. 


MSp  and  Particles.  All  samples  from  Cuitzeo  contain  a  variable  mix 
of  MSp,  magnetic  grains,  framboidal  spherules,  and/or  weakly 
magnetic  volcanic  glass  (Fig.  4;  SI  Appendix,  Tables  4  and  5).  MSp 
range  from  25-100  pm  in  diameter,  averaging  60  pm  (Fig.  5),  and 
typically  appear  as  highly  reflective,  black  spheroids  (Fig.  5  A 
and  B;  SI  Appendix,  Fig.  6  B  and  E),  although  shapes  such  as 
ovals,  doublets,  dumbbells,  and  tear-drops  frequently  occur  (Fig.  5 
D-G;  SI  Appendix,  Fig.  6^4,  C,  andD).  The  MSp  are  conspicuous 
and  abundant  at  2.8  m  (2,000/kg),  where  they  form  a  sharp  peak 
in  the  YDB.  They  average  approximately  100  MSp /kg  above  the 
YDB,  and  none  were  detected  at  depths  of  3.6  to  3.0  m. 

SEM  observation  of  outer  surfaces  of  all  MSp  analyzed  reveals 
a  surficial  crystalline  pattern  that  is  dendritic  (Fig.  5  D  and  F)  or 
polygonal  like  a  soccer-ball  (Fig.  5/4;  SI  Appendix,  Fig.  6 E).  These 
patterns  are  indicative  of  melting  with  rapid  quenching11,  which 
precludes  diagenetic,  biogenic,  or  detrital  origins.  Several  MSp 
clearly  display  evidence  of  interspherule  collisions  while  solidify¬ 
ing,  causing  fusion  (Fig.  5  D  and  E ).  EDS  analyses  of  the  recov¬ 
ered  MSp  show  that  they  are  comprised  of  magnetite  (Fe  oxide 
>96%)  with  very  low  abundances  of  other  elements  (SI  Appendix, 
Table  4).  Because  titanium  is  present  at  only  trace  levels  in  the 
MSp,  they  are  not  titanomagnetite  grains,  which  are  ubiquitous 
throughout  the  Cuitzeo  core. 

We  observed  variable  abundances  of  irregularly  shaped  titano¬ 
magnetite  grains  throughout  the  sequence,  peaking  well  below 
the  YD  onset  (Fig.  4).  Typically  black  and  moderately  reflective, 
these  grains  are  monocrystalline  and  geochemically  distinct  from 
MSp  and  are  interpreted  to  be  detrital  material  from  the  local 
volcanic  regolith.  Pinter  et  al.  (27)  speculated  that  all  previously 
reported  YDB  MSp  are  actually  detrital  grains  or  framboidal 
spherules,  but  this  is  refuted  at  Lake  Cuitzeo  because  such  grains 
are  readily  distinguishable  from  YDB  MSp  by  the  quench-melted 
textures  visible  by  SEM  imaging. 

Irregularly  shaped  volcanic  glass  grains  (<  1  pm  to  several  mm 
long)  occur  in  all  layers  and  contain  up  to  10%  Fe  causing  them  to 
be  magnetic.  They  were  observed  in  relatively  uniform  abun¬ 
dances  throughout  the  section,  uncorrelated  with  the  MSp  peak. 
Most  appear  dark-gray  to  black,  translucent-to-opaque,  and  are 
highly  reflective.  SEM  imaging  revealed  vesicles  (gas  bubbles) 
that  formed  in  the  glass  as  it  cooled.  Using  a  standard  igneous 


nPetaev  Ml,  Jacobsen  SB,  Basu  AR,  Becker  L,  Magnetic  Fe,  Si,  Al-rich  impact  spherules 
from  the  P-T  boundary  layer  at  Graphite  Peak,  Antarctica,  35th  Lunar  and  Planetary 
Conference,  March  15-19,  2004,  Houston,  TX,  1216  (abstr.). 
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rock  binary  diagram,  the  glass  grains  are  identified  as  alkali 
basalt,  andesite,  rhyolite,  rhyodacite,  tephrite,  and  latite  (tra- 
chyandesite)  (SI Appendix,  Fig.  7).  These  glass  grains  are  consid¬ 
ered  to  be  of  local  volcanic  origin.  In  contrast,  none  of  the 
quenched  MSp  plotted  on  the  graph,  indicating  a  different  origin. 

Framboidal  spherules  (approximately  20  to  60  pm,  averaging 
approximately  40  pm  in  diameter)  are  concentrations  of  euhedral 
microcrystals  typical  of  diagenetic  sedimentary  pyrite  (28) 
(Fig.  5C;  SI  Appendix,  Fig.  6 F).  Framboidal  spherules  were  found 
intermixed  with  MSp,  forming  a  distinct  abundance  peak  in  the 
YDB  layer  at  2.75  m,  where  they  occur  at  10,000/kg  (Fig.  4). 
None  were  detected  below  the  base  of  the  YDB  at  2.80  m.  On 
the  other  hand,  above  2.75  m,  framboidal  spherules  were  rela¬ 
tively  common,  averaging  a  few  hundred  per  kg.  The  origin  of 
the  framboidal  spherules  is  unclear,  but  their  formation  most 
likely  resulted  from  the  onset  of  anoxic  conditions  beginning  near 
12.9  ka. 

We  have  compared  and  contrasted  the  geochemistry  of  MSp, 
framboidal  spherules,  glassy  grains,  and  magnetic  grains  from  the 
Cuitzeo  YDB  layer  using  SEM-EDS.  A  ternary  diagram  in  Fig.  6A 
compares  FeO,  Si02,  and  Ti02  and  demonstrates  that  the  MSp 
are  geochemically  dissimilar  to  volcanogenic  material.  Next,  the 
geochemistry  of  the  YDB  MSp  was  plotted  on  ternary  diagrams 
and  compared  to  that  of  other  types  of  spherules  and  melted  ma¬ 
terial.  In  Fig  6 B,  no  similarity  was  observed  when  compared  with 
particles  representative  of  cosmic  influx  and  meteoritic  ablation, 
including  >700  meteorites  and  cosmic  spherules  (SI  Appendix, 
Table  6).  Instead,  geochemical  values  for  Cuitzeo  MSp  are  similar 
to  those  for  >1,000  tektites  (glassy,  melted  impact  material)  and 
MSp  from  11  craters/strewnfields  formed  by  ET  impact  into 
terrestrial  rocks  (Fig  6C;  SI  Appendix,  Table  6),  suggesting  that 
Cuitzeo  MSp  also  formed  by  cosmic  impact. 

NDs.  A  definitive  proxy  for  the  YD  impact  event  is  an  abundant 
and  diverse  assemblage  of  NDs  within  the  YDB  layer  across 
North  America  (2,  3)  and  in  the  Greenland  Ice  Sheet  (4).  Four 
allotropes  of  NDs  have  been  previously  identified,  including 
cubic  NDs  and  lonsdaleite,  which  has  only  been  found  on  Earth 
associated  with  impact  craters  or  within  meteoritic  material  (29- 
31). The  other  two  allotropes,  n-diamond  and  i-carbon,  were  first 
discovered  in  the  laboratory  and  are  considered  to  be  either  dia- 
mond-like  carbon  or  modified  cubic  NDs  (32,  33). 

In  a  comprehensive  investigation  of  NDs  at  Cuitzeo,  we  se¬ 
lected  eight  bulk  samples  spanning  1.4  m  from  3.6  to  2.2  m  that 
included  the  YD  interval  within  an  age  range  from  21  to  10  ka. 
Noncontinuous  approximate  1-cm-thick  samples,  ranging  from 
approximately  4  to  23  g  (approximately  3  to  20  cm3)  were  taken 
at  intervals  of  5  cm  across  the  inferred  YDB  and  at  intervals  of  10 
to  50  cm  away  from  the  YDB.  Using  procedures  described  in 
Kennett  et  al.  (2,  3),  we  produced  acid-resistant  residue  for  each 
bulk  sample  to  concentrate  any  NDs  (Methods',  SI  Appendix, 
Methods).  In  response  to  a  challenge  by  Daulton  et  al.  (8)  that 


the  YDB  contains  no  NDs,  we  examined  these  residues  using 
a  highly  comprehensive  suite  of  analyses  surpassing  previous 
investigations.  Analyses  were  conducted  on  high-resolution  trans¬ 
mission  electron  microscopes  (HRTEM)  and  scanning  trans¬ 
mission  electron  microscopes  (STEM),  discussed  below. 

These  analyses  clearly  identified  a  single  major  peak  in  NDs, 
centered  across  two  samples  at  2.8  and  2.75  m  within  the  dark 
layer  dating  to  the  YD  onset.  NDs  were  rare  below  the  2.9  m  layer 
[<1  parts  per  billion  (ppb)],  whereas  above  the  2.75  m  layer, 
NDs  were  observed  at  low  levels  from  4  to  10  ppb,  likely  due  to 
reworking.  Initial  examination  of  the  2.8  m  sample  using  STEM 
and  HRTEM  revealed  a  striking  panorama  of  tens  of  thousands 
of  nanocrystalline  carbon  particles  ranging  in  shape  from  sphe¬ 
rical  to  elongate  to  euhedral  (SI  Appendix,  Fig.  8).  These  particles 
varied  in  diameter  from  approximately  1  to  10  nm,  averaging  ap¬ 
proximately  4  nm,  and  were  typically  embedded  in  amorphous 
carbon,  as  Tian  et  al.  (14)  described.  We  identified  three  of  four 
previously  reported  ND  variants,  of  which  n-diamond  was  most 
abundant,  with  lesser  amounts  of  i-carbon  and  lonsdaleite.  The 
presence  and  concentrations  of  cubic  NDs  are  unclear,  for  rea¬ 
sons  discussed  below.  The  NDs  exhibit  maximum  abundance  at 
2.8  m  of  approximately  100  ±  50  ppb  (SI  Appendix,  Table  5),  si¬ 
milar  to  the  estimates  for  YDB  NDs  across  North  America  and 
Greenland  (2-4). 

EDS  analyses  of  five  samples  indicate  that  the  observed  crys¬ 
tals  are  dominantly  carbon  (averaging  96  atomic  %)  with  minor 
concentrations  of  other  elements  (SI  Appendix,  Table  7).  Using 
noncopper  grids  (gold  and  molybdenum)  for  all  five  EDS  ana¬ 
lyses,  we  detected  no  copper,  refuting  speculation  by  Daulton 
et  al.  (8)  that  YDB  NDs  are  comprised  of  copper  and  not  carbon. 

Selected  area  diffraction  (SAD)  was  used  to  identify  the  domi¬ 
nant  types  of  carbon  crystals  present  in  the  samples.  Many  ana¬ 
lyses  revealed  ring  diffraction  patterns  with  d-spacings  that  are 
consistent  with  n-diamond  (Fig.  1A;  SI  Appendix,  Table  8),  with 
a  face-centered  cubic  structure  with  space  group  Fm3m  (34). 
Other  SADs  indicated  d-spacings  characteristic  of  i-carbon 
(35)  (Fig.  1B\  SI  Appendix,  Table  8),  a  primitive  cubic  crystalline 
structure  with  space  group  P2X3  or  F4232  (32).  Crystals  with  d- 
spacings  consistent  with  cubic  nanodiamond  were  also  observed, 
but  identification  was  not  conclusive. 

Fast  Fourier  transform  (FFT)  is  an  analytical  procedure  that 
produces  diffraction  patterns  for  small  single  crystals.  HRTEM 
and  FFT  images  of  the  nanocrystals  reveal  lattices  and  d-spacings 
that  are  consistent  with  lonsdaleite,  space  group  F^/mmc 
(Fig.  8  A  and  B;  SI  Appendix,  Table  8)  and  with  n-diamond,  space 
group  Fm3m  (Fig.  8  C  and  D). 

Electron  energy  loss  spectroscopy  (EELS)  was  utilized  to  con¬ 
firm  that  the  observed  nanoparticles  are  carbon  and  crystalline,  as 
previously  indicated  by  HRTEM,  FFT,  and  SAD  analyses.  EELS 
spectra  further  indicate  that  the  nanoparticles  display  sp3  and  sp2 
bonding  consistent  with  n-diamonds  and  i-carbon.  Independent 
published  spectra  of  NDs  are  plotted  for  comparison  with  those 


Fig.  6.  Ternary  geochemical  diagrams:  (A)  Cuitzeo  magnetic  impact  spherules  compared  to  volcanogenic  titanomagnetite  and  glassy  grains,  as  well  as  fram¬ 
boidal  spherules.  Cuitzeo  magnetic  impact  spherules  are  nonvolcanogenic.  Of  two  framboidal  spherules  analyzed,  one  overlaps  the  magnetic  impact  spherules 
and  one  does  not.  Neither  exhibits  quench  melting.  ( B )  Cosmic  particles  compared  to  Lake  Cuitzeo  magnetic  impact  spherules,  indicating  they  are  noncosmic. 
(C)  Terrestrial  impact  materials  compared  to  Lake  Cuitzeo  magnetic  impact  spherules,  showing  close  geochemical  match. 
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Fig.  7.  SAD  patterns  of  NDs  from  2.7  m.  (A)  D-spacings  indicative  of 
n-diamonds.  ( B )  D-spacings  indicative  of  i-carbon. 

for  Cuitzeo  NDs  in  Fig.  9  (4,  32,  34,  36,  37).  The  small  prepeak  at 
approximately  284  eV  indicates  sp2  bonding  and  is  similar  for  all 
curves.  This  prepeak  has  been  previously  interpreted  by  others  to 
result  from  a  graphitic  or  amorphous  coating  on  the  NDs  (SI 
Appendix ,  EELS  Analyses).  The  spectrum  does  not  permit  defi¬ 
nitive  identification  of  the  other  two  allotropes,  cubic  NDs, 
and  lonsdaleite.  Another  EELS  spectrum  confirmed  a  lack  of 
copper  in  the  analyzed  particles  (SI  Appendix,  Fig.  9). 

Energy-filtered  TEM  analysis  (EFTEM),  an  EELS-related 
mapping  technique,  was  used  to  investigate  the  elemental  com¬ 
position  of  the  nanocrystals.  We  produced  two  maps  of  a  region 
of  nanoparticles  by  searching  for  specific  energy  levels  indicative 
of  carbon.  First,  a  “zero-loss”  image  (Fig.  1CL4)  was  acquired  and 
inverted  to  negative  for  easier  comparison.  It  displays  brighter 
particles  (numbers  1  through  4)  that  were  embedded  in  amor¬ 
phous  carbon,  but  exhibit  a  crystalline  structure.  Second,  a  “jump 
ratio”  image  (Fig.  10 B)  shows  the  same  particles  (numbers  1 
through  4).  Their  brighter  shading  indicates  the  presence  of 
sp3  bonding  characteristic  of  crystalline  carbon,  including  NDs. 
These  results,  coupled  with  the  previous  SAD  and  FFT  analyses, 
are  consistent  with  NDs  and  inconsistent  graphite,  graphene,  and 
copper  (see  SI  Appendix,  EFTEM  Analysis). 


Fig.  8.  Crystallographic  data  for  NDs  from  2.8  m.  (A)  HRTEM  image  of  mono¬ 
crystalline  nanoparticle  identified  as  lonsdaleite.  The  (101)  and  (110)  planes 
are  visible  with  d-spacings  of  1.93  and  2.18  A),  respectively.  ( B )  FFT  of  same 
lonsdaleite  crystal  above.  The  values  adjacent  to  each  spot  indicate  the  re¬ 
ciprocal  lattice  vector.  Image  reveals  (101)  planes  with  lattice  spacing  of 
1 .93  A,  consistent  with  lonsdaleite.  (O  HRTEM  image  displaying  typical  lattice 
spacing  of  n-diamond.  The  1 .78  A  measurement  represents  the  (200)  planes, 
consistent  with  n-diamond.  (D)  FFT  of  same  n-diamond  shown  above. 
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Fig.  9.  EELS  spectra  for  NDs  from  2.8  m.  A  typical  carbon  peak  of  approxi¬ 
mately  295-300  eV  shows  that  the  particle  is  carbon.  Published  n-diamond 
and  i-carbon  spectra  (dotted  lines)  are  shown  for  comparison  (4,  32,  36,  37). 
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Twinned  NDs.  These  NDs  are  made  up  of  two  or  more  crystals  that 
share  a  common  lattice  plane  (the  twin  plane)  and  grow  symme¬ 
trically  away  in  different  orientations  (Fig.  11).  Twinning  is  com¬ 
monly  observed  in  commercial  NDs  formed  by  carbon  vapor 
deposition  (CVD),  during  which  NDs  crystallize  from  gaseous 
carbon,  typically  at  high  temperatures  in  an  inert  atmosphere. 
Twinned  cubic  NDs  are  common  in  meteorites  (38),  having 
formed  in  space  through  a  process  possibly  analogous  to  CVD 
(39).  They  also  are  found  in  impact  craters  (39)*$,  where  they 
formed  upon  impact  from  terrestrial  carbon.  Twinned  lonsdaleite 
has  been  observed  in  meteorites  and  associated  with  impact  cra- 
terstt.  Twins  can  form  in  numerous  configurations,  including  “ac¬ 
cordion  twins,”  which  exhibit  folded,  pleat-like  lattice  planes,  and 
fivefold  “star  twins”  (38),  as  observed  by  Tian  et  al.  (14)  in  the 
YDB  layer  from  Lommel,  Belgium.  At  Cuitzeo,  most  NDs  were 
twinned  n-diamond  and  i-carbon  and  only  occasionally  were 
monocrystalline  NDs  observed.  Twinned  lonsdaleite  with  d-spa- 
cings  of  2.06  A  and  1.93  A  was  observed  occasionally  (Fig.  11 B). 

Cubic  NDs.  Cubic  NDs  were  previously  identified  in  the  YDB 
(2-4),  and  subsequently,  Tian  et  al.  (14)  confirmed  cubic  NDs 
in  the  YDB.  In  the  Cuitzeo  section,  however,  we  could  not  un¬ 
equivocally  identify  the  cubic  allotrope.  This  may  be  due  to  mask¬ 
ing  by  i-carbon  and/or  n-diamonds,  which  share  some  d-spacings 
with  cubic  NDs  (SI  Appendix,  Table  8).  Also,  cubic  NDs  possess 
so-called  “forbidden  reflections,”  such  as  the  1.78  A  d-spacing, 
that  are  typically  invisible  in  cubic  SAD  patterns  but  are  some¬ 
times  apparent  in  twinned  cubic  NDs,  most  likely  due  to  double 
diffraction  (38).  Because  n-diamonds  also  display  these  forbidden 
reflections,  twinned  n-diamonds  cannot  be  easily  differentiated 
from  twinned  cubic  NDs.  Thus,  it  is  possible  that  some  of  the  ap¬ 
parent  n-diamonds  from  Cuitzeo  are  actually  twinned  cubic  NDs. 

Carbon  Onions  and  Ribbons.  In  the  YDB  at  Lommel,  Belgium,  Tian 
et  al.  (14)  observed  carbon  onions,  nanonmeter-sized  nanoparti¬ 
cles  constructed  of  concentric  carbon  shells.  Carbon  onions 
typically  are  comprised  of  three  to  eight  closed  shells,  ranging 
from  3  to  100  nm  wide  (40),  and  may  be  either  hollow  or  contain 
NDs.  Broken  and  whole  carbon  onions,  ranging  from  ovoid  to 
round  and  from  2  to  10  nm  wide,  were  observed  in  the  2.75  m 
sample  using  HRTEM  imaging  (Fig.  12;  SI  Appendix,  Fig.  10). 
These  display  carbon  shells  with  variable  spacings,  ranging  from 
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Fig.  10.  EFTEM  maps  of  NDs.  (A)  An  inverted  "zero-loss"  image  displaying 
brighter  nanoparticles  (numbered).  Inset  shows  4-nm-wide  nanoparticle 
number  1  that  exhibits  crystalline  lattice  spacings  of  1.30  and  2.12  A,  consis¬ 
tent  with  i-carbon.  Resolution  is  low  due  to  surrounding  amorphous  carbon. 
(B)  A  map  of  particles  detected  with  characteristic  carbon  signature  at 
299  ±  5  eV.  Brighter  particles  at  numbers  1  through  4  correspond  to  the  Left 
panel,  indicating  sp3  bonding  typical  of  NDs.  Darker  area  near  number  3  is 
due  to  sp2  bonding  in  the  amorphous  carbon  film. 

approximately  3.47  to  4.37  A  with  an  average  of  4.07  A.  Although 
it  is  not  possible  to  identify  the  exact  carbon  allotrope,  these  spa¬ 
cings  are  somewhat  consistent  with  fullerenes  (4.09  A),  which 
were  previously  detected  in  the  YDB  layer  in  North  America 
(1).  About  half  of  the  carbon  onions  enclose  apparent  polycrystal¬ 
line  i-carbon,  as  measured  by  FFT.  However,  definitive  analysis  is 
not  possible  because  the  NDs  are  very  small  and  surrounded  by 
amorphous  carbon  and  carbon  onion  shells. 

Because  natural  wildfires  can  produce  carbon  onions,  Tian  et 
al.  (14)  postulated  that  YDB  NDs  may  have  formed  during  wild¬ 
fires  by  high-pressure  compression  within  carbon  onions.  How¬ 
ever,  there  is  no  evidence  to  support  that  suggestion  because 
no  naturally  formed  carbon  onions  have  ever  been  found  to  con¬ 
tain  NDs  (SI Appendix,  Fig.  11  and  Carbon  Onions).  On  the  other 
hand,  diamond-bearing  carbon  onions  are  known  to  form  during 
production  of  NDs  by  trinitrotoluene  (TNT)  detonation  (41)  un¬ 
der  high  temperatures  and  oxygen-deficient  conditions  that  are 
known  to  occur  during  an  impact  event,  as  discussed  below. 

We  also  observed  carbon  ribbons,  which  are  multinanometer- 
long  parallel  chains  of  carbon  atoms  found  in  interstellar  dust 
particles,  or  IDPs,  and  meteorites  (42)  (SI  Appendix ,  Fig.  10). 
They  have  also  been  observed  in  nanodiamond-rich  residues  from 
TNT  detonation  (41). 


Fig.  11.  HRTEM  images  of  twinned  NDs  from  the  2.8  m  layer.  Double  yellow 
lines  represent  lattice  planes  and  the  numbers  indicate  d-spacings  in  A. 
Arrows  are  parallel  to  common  twinning  plane.  (A)  Star-twin  ND  with  five¬ 
fold  star-like  morphology.  ( B )  Accordion  twin  lonsdaleite  with  pleated  mor¬ 
phology.  (C)  Twin  with  multiple  folds.  (D)  "Scalloped"  twin. 


Fig.  12.  Carbon  onions  from  2.75  m.  (A)  HRTEM  image  displays  10  nm  car¬ 
bon  onion.  Parallel  strands  nearby  appear  to  be  carbon  ribbons.  (B)  Drawing 
illustrates  shells  of  carbon  onion  and  nanoparticle  shown  in  A.  (Q  FFT  of 
enclosed  crystal  with  d-spacings  generally  consistent  with  i-carbon  but  with 
insufficient  resolution  to  be  definitive. 

Discussion 

Multiple  hypotheses  have  been  proposed  to  explain  YDB  mar¬ 
kers.  Nearly  all  can  be  rejected,  as  follows: 

"Cosmic  Rain."  Several  workers  have  speculated  that  YDB  MSp 
and  NDs  represent  the  normal  “cosmic  rain”  that  falls  continu¬ 
ously  from  space  (6,  13,  27).  This  is  a  testable  hypothesis,  but 
those  authors  conducted  no  such  tests.  For  Cuitzeo,  we  analyzed 
all  available  MSp  geochemically  and  compared  them  to  >700  me¬ 
teorites  and  cosmic  spherules  from  Greenland,  Antarctica,  and 
elsewhere  (Fig.  6).  The  MSp  geochemistry  indicates  that  they 
are  noncosmic  and,  instead,  appear  to  be  melted  terrestrial  ma¬ 
terial.  Furthermore,  the  high  abundance  of  YDB  MSp  indicates 
that  cosmic  rain  is  an  unlikely  cause  because  cosmic  spherules  in 
sediment  are  very  rare.  For  example,  the  1991  EUROMET  Ant¬ 
arctic  micrometeorite  collection  from  continuously  deposited  ice 
contains  only  0.015  microspherules  (50-100  pm)  per  kg  of  ice 
(1  L)  (43).  In  contrast,  Cuitzeo  contains  >2,000  MSp  per  kg,  to¬ 
talling  130,000  x  more  for  the  same  sample  weight. 

Cosmic  NDs  occur  in  meteorites  and  cosmic  dust,  but  Tian 
et  al.  (14)  concluded  that  YDB  NDs  are  not  cosmic  because  they 
display  <513C  abundances  (—28.1  to  —  26.3%c)  that  are  terrestrial. 
Their  work  confirmed  the  results  of  Bunch  et  al.§§,  who  reported 
isotopic  analyses  on  diamond-rich  residue  from  12  YDB 
sites  that  exhibited  terrestrial  values  for  <513C  (—30  to  —22 %o) 
and  <515N  (—17  to  15%c).  Furthermore,  these  values  are  similar 
to  those  for  NDs  in  the  Cretaceous-Paleogene  boundary  layer 
(KPg,  formerly  K-T),  where  Gilmour  et  al.  (44)  measured  <513C 
(—18 %o)  and  <515N  (6%o)  and  concluded  that  KPg  NDs  were 
formed  during  impact  from  terrestrial  material,  most  likely  from 
carbon-rich  target  rocks.  Isotopic  analyses  of  the  carbon-rich 
YDB  interval  at  Cuitzeo  yielded  values  ranging  from  —23  to 
—19 %c  for  <513C  consistent  with  the  formation  of  Cuitzeo  NDs 
from  terrestrial,  not  cosmic,  carbon. 

Wildfires.  Typical  wildfire  temperatures  range  from  approximately 
900  to  1,200  °C  with  the  highest  reported  being  1,450  °C  (45). 
These  temperatures  are  too  low  to  melt  magnetite  into  MSp 
(1,540  °C)  and  too  high  for  the  survival  of  NDs,  which  combust 
between  400  and  600  °C  (39).  Neither  proxy  has  been  reported  in 
normal  wildfires.  This  conclusion  is  reinforced  by  our  observation 
that  the  Cuitzeo  3.1  m  peak  in  charcoal  is  indicative  of  a  major 
wildfire  episode,  and  yet,  displays  low  levels  of  MSp,  CSp,  or 
NDs  (Fig.  4;  SI  Appendix ,  Table  5).  Another  indication  that 
the  YDB  proxies  are  not  wildfire-related  is  that  marker  peaks 
(2.80  to  2.75  m)  were  deposited  earlier  than  the  wildfire  charcoal 
peak  (2.70  to  2.65  m). 
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Volcanism.  Low-energy  volcanism  produces  silicate  spherules  but 
does  not  distribute  them  widely  (46).  Conversely,  high-energy 
eruptions  capable  of  disseminating  ejecta  widely  do  not  appear 
to  produce  spherules.  For  example,  the  Toba  eruption  at  approxi¬ 
mately  75  ka,  one  of  the  largest  of  the  last  5  million  years,  ejected 
debris  up  to  2,000  km,  yet  no  spherules  have  been  detected  (46). 
Likewise,  we  analyzed  tephra  samples  from  the  Laacher  See  erup¬ 
tion  layer  in  Germany  that  preceded  the  YD  onset  and  found  no 
MSp  or  NDs.  On  the  other  hand,  the  Cuitzeo  YDB  layer  contains 
almost  no  volcanic  material;  the  closest  volcanic  layers  are  approxi¬ 
mately  18  kyr  earlier  and  approximately  8  kyr  later  (15),  indicating 
that  no  significant  local  volcanic  eruptions  occurred  near  12.9  ka. 

Anthropogenesis.  NDs  have  never  been  reported  in  industrial  by¬ 
products,  and  because  MSp,  CSp,  and  NDs  are  deeply  buried  at 
approximately  3  m  in  well-stratified  deposits  at  Cuitzeo,  there  was 
no  reasonable  chance  for  human  contamination. 

Potential  Misidentification  of  Markers.  Surovell  et  al.  (6)  reported 
finding  no  YDB  MSp  peaks,  although  claiming  to  follow  the  pro¬ 
tocol  of  Firestone  et  al.  (1)  for  quantification  of  MSp,  and  con¬ 
cluded  that  Firestone  et  al.  misidentified  and/or  miscounted  the 
MSp.  Later,  Lecompte  et  al.1™  independently  examined  two  YDB 
sites  common  to  Firestone  et  al.  and  Surovell  et  al.  They  reported 
that  “spherule  abundances  are  consistent  with  those  of  Firestone 
et  al.”  and  “inconsistent  with  the  results  of  Surovell  et  al.”  They 
also  concluded  that  Surovell  et  al.  altered  the  prescribed  MSp 
protocol  in  fatal  ways,  particularly  by  not  observing  requirements 
for  sample  thickness,  sample  weight,  and  size  sorting.  We  consid¬ 
er  these  discrepancies  significant  enough  to  negate  the  conclu¬ 
sions  of  Surovell  et  al.  (SI Appendix,  Surovell  et  al). 

Daulton  et  al.  (8)  found  no  YDB  NDs  at  Arlington  Canyon, 
California,  or  at  Murray  Springs,  Arizona,  as  earlier  reported  in 
Kennett  et  al.  (2,  3).  They  searched  for  NDs  in  “microcharcoal 
aggregates”  from  the  Murray  Springs  YDB  site  and,  finding  none, 
claimed  to  refute  the  previous  results.  However,  Kennett  et  al. 
never  claimed  to  find  NDs  in  charcoal,  and  instead,  observed 
NDs  at  Murray  Springs  in  acid-resistant  residues  from  bulk  sedi¬ 
ment  (2,  3),  which  Daulton  et  al.  did  not  investigate. 

Daulton  et  al.  (8)  further  speculated  that  Kennett  et  al.  (2,  3) 
misidentified  YDB  NDs,  observing  copper  instead,  which  dis¬ 
plays  d-spacings  nearly  identical  to  n-diamond  and  i-carbon.  In 
addition,  Daulton  et  al.  pointed  out  that  graphene  and/or  gra- 
phane  have  d-spacings  similar  to  lonsdaleite  and  that  the  lonsda¬ 
leite  diffraction  pattern  reported  from  Arlington  Canyon  by  the 
Kennett  et  al.  (2)  was  missing  the  lonsdaleite  diffraction  line  at 
1.93  A.  However,  in  YD-aged  ice  in  Greenland,  Kurbatov  et  al. 
(4)  identified  lonsdaleite  with  the  1.93-  A  line,  which  definitively 
demonstrates  that  those  Greenland  nanoparticles  cannot  be  gra¬ 
phene  or  graphane.  At  Lake  Cuitzeo,  numerous  NDs  have  been 
identified  with  the  1.93  A  (101)  line,  as  shown  in  Fig.  8^4and  B 
and  Fig.  11 B,  eliminating  the  possibility  that  these  crystals  are 
graphene  or  graphane.  SAD  and  all  other  analyses  conclusively 
show  that  the  Cuitzeo  nanoparticles  analyzed  have  d-spacings 
consistent  with  lonsdaleite  and  other  NDs.  In  independent  sup¬ 
port  of  NDs  in  the  YDB,  Tian  et  al.  (14)  and  Van  HoeseF  iden¬ 
tified  cubic  NDs  in  the  YDB  layer  in  Europe. 

Regarding  CSp,  Scott  et  al.  (7)  speculated  that  those  found  at 
YDB  sites  (1-3)  are  simply  charred  fungal  sclerotia,  which  are 
ball-like  clusters  of  long,  branching  filamentous  structures,  com¬ 
mon  to  some  fungi.  The  CSp  from  Cuitzeo  and  other  YDB  sites 
are  unmistakably  different  from  sclerotia  in  numerous  critical 
characteristics.  In  particular,  charred  and  uncharred  sclerotia 
have  textured,  filamentous,  low-reflectivity  interiors,  whereas 
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at  Cuitzeo,  SEM  imaging  demonstrates  that  CSp  have  smooth, 
glassy,  highly  reflective  interiors  with  no  evidence  of  filamentous 
structure  observed  in  fungal  sclerotia  (or  cellular  structure  found 
in  charcoal)  (SI  Appendix,  Fig.  5). 

Cuitzeo  CSp  also  contain  numerous  noncarbon  particles,  in¬ 
cluding  aluminosilicates,  indicating  that  these  cannot  be  primary 
biological  entities,  such  as  sclerotia.  In  support  of  this,  several 
lines  of  evidence  support  the  formation  of  CSp  during  biomass 
burning.  For  example,  Firestone  et  al.  (1)  reported  the  produc¬ 
tion  of  CSp  in  modern  wildfires,  and  laboratory  experiments  have 
demonstrated  the  production  of  CSp  from  charred  tree  resin  at 
approximately  500  °C***.  These  CSp  are  morphologically  iden¬ 
tical  to  those  found  in  the  YDB  but  contain  no  NDs.  Also,  CSp 
similar  to  those  found  in  the  YDB  have  been  reported  by  Harvey 
et  al.  (47),  who  observed  vesicular  CSp  in  the  impact  layer  at  the 
KPg,  and  suggested  that  CSp,  along  with  aciniform  soot,  formed 
during  impact  into  carbon-rich  target  rocks. 

Cosmic  Impact  as  Only  Viable  Hypothesis 
Impact-related  CVD.  Tian  et  al.  (14)  speculated  that  YDB  NDs 
formed  by  CVD,  although  they  offered  no  details.  In  the  labora¬ 
tory,  formation  of  NDs  by  CVD  requires  intense  heating  of  car¬ 
bon  vapor  within  an  inert  atmosphere,  conditions  not  known  to 
exist  naturally  at  Earth’s  surface  [SI  Appendix,  CVD  ( Carbon 
Vapor  Deposition)].  ET  impacts  are  the  only  known  natural  events 
capable  of  generating  CVD -like  conditions  under  a  reduced-oxy¬ 
gen  atmosphere  (39).  This  CVD  mechanism  has  been  proposed 
for  the  KPg,  where  <513C  and  <5 15  N  values  for  cubic  NDs  suggest 
they  formed  from  carbon  that  is  terrestrial  and  not  cosmic  (44). 

Comets.  Based  upon  astrophysical  observations  and  modeling, 
Napier  (48)  proposed  that  YDB  impact  markers  were  produced 
when  Earth  encountered  a  dense  trail  of  material  from  a  large 
already  fragmented  comet.  His  model  predicts  cluster  airbursts 
and/or  small  cratering  impacts  that  could  account  for  the  wide 
distribution  of  YD  impact  debris  across  more  than  10%  of  the 
planet,  including  Cuitzeo.  Most  comets  eventually  break  up  as 
they  transit  the  inner  solar  system,  and  previously  unknown  frag¬ 
mented  comets  are  discovered  by  space-borne  telescopes,  such  as 
the  Solar  and  Heliospheric  Observatory,  on  average  every  4  y.  As 
evidence,  Earth  is  bombarded  at  an  average  rate  of  once  every  5  d 
by  one  of  72  meteor  streams  or  “showers,”  massive  clouds  of  deb¬ 
ris  from  fragmented  comets.  These  well-known  meteor  showers, 
e.g.,  Perseids,  Geminids,  Taurids,  etc.,  are  highly  dispersed,  but  in 
the  recent  geologic  past,  each  stream  was  far  more  condensed, 
containing  many  large,  potentially  destructive  fragments.  Cur¬ 
rently,  the  Taurid  Complex  contains  19  large  near-earth  Apollo 
asteroids,  with  diameters  ranging  from  approximately  1.5  km 
(6063  Jason)  to  approximately  5  km  (4184  Cuno)  (48).  None 
of  these  currently  threatens  Earth  but  may  do  so  in  the  future. 

Impact  Dynamics.  Earth  has  been  subjected  to  a  continuous, 
although  intermittent  bombardment  by  impactors  with  diameters 
ranging  from  microns  to  tens  of  kilometers;  velocities  range  from 
approximately  1 1  km/s  to  73  km/s  with  typical  values  of  17  km/s 
for  asteroids  and  51  km/s  for  comets.  The  term  “cosmic  impact” 
evokes  images  of  craters  ranging  from  the  50-kyr-old,  1-km-dia- 
meter  Meteor  Crater  to  the  2-billion-year-old,  200-km-diameter 
Vredefort  crater  (49,  50).  For  these  crater-forming  events  that 
have  peak  impact  pressures  in  the  range  of  hundreds  of  GPa, 
impact  dynamics  and  shock  wave  metamorphic  effects  are  well 
understood  (49,  50).  An  ET  impact  is  the  only  natural  mechanism 
known  to  produce  major  coeval  abundances  in  cubic  NDs,  lons- 


***Kimbel  D,  West  A,  Kennett  JP,  A  new  method  for  producing  nanodiamonds  based  on 
research  into  the  Younger  Dryas  extraterrestrial  impact,  AGU  Fall  Meeting,  December 
15-19,  2008,  Eos  Trans.  AGU,  90(52),  Fall  Meet.  Suppl.,  #PP13C-1470. 
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daleite,  and  quench-melted  MSp,  both  of  which  co-occur  in  im¬ 
pact  events,  including  Ries  crater  and  the  KPg  (39). 

Based  on  hundreds  of  shock-recovery  experiments  by  one  of 
the  authors  of  this  article  (DeCarli),  the  formation  of  lonsdaleite 
in  graphite-bearing  gneisses  in  the  Ries,  Popigai,  and  other  im¬ 
pact  craters  is  in  complete  accord  with  static  high-pressure  data 
on  solid-solid  transformation  of  graphite  to  lonsdaleite  and  cubic 
NDs  (29-31).  However,  this  transformation  does  not  readily  ex¬ 
plain  the  NDs  found  at  the  KPg  boundary  or  in  the  YDB.  Based 
on  available  evidence,  it  seems  unlikely  that  the  YDB  NDs 
formed  by  shock  compression  of  terrestrial  graphite,  and  instead, 
our  preferred  mechanism  invokes  the  interaction  of  an  ET object 
with  Earth’s  atmosphere  (49).  If  incoming  objects  are  relatively 
small,  virtually  all  kinetic  energy  is  transferred  to  the  atmosphere 
at  high  altitudes,  creating  an  air  shock  with  temperatures  up  to 
tens  of  thousands  degrees  Kelvin.  These  are  the  familiar  shooting 
stars,  the  remains  of  which  may  be  collected  as  cosmic  dust. 
Although  shock  pressures  due  to  solid-air  interaction  are  modest 
at  high  altitudes,  larger  objects  may  be  disrupted  and  fragmented 
as  pressure  builds  due  to  increasing  atmospheric  density  at  lower 
altitudes.  This  breakup  is  especially  likely  if  the  object  was  loosely 
consolidated  or  low  density  like  a  comet.  When  an  incoming  ET 
object  encounters  the  atmosphere  and  breaks  apart,  individual 
pieces  rapidly  decelerate  due  to  the  marked  increase  in  the  ratio 
of  cross-sectional  area  to  mass.  Area  of  the  luminous  air  shock  is 
correspondingly  increased,  with  the  result  that  the  object  appears 
to  “explode”  in  a  fireball.  For  an  object  traveling  at  30  km/s,  air 
shock  pressure  would  be  approximately  20  MPa  at  20  km  altitude, 
approximately  170  MPa  at  10  km,  and  approximately  900  MPa  at 
sea  level.  For  an  air  shock  of  170  MPa,  the  pressure  exceeds  un¬ 
confined  compressive  strengths  of  many  rocks. 

These  energetic  events  are  often  termed  “atmospheric  im¬ 
pacts”  to  distinguish  them  from  more  familiar  crater-forming 
events.  For  example,  the  craterless  Tunguska  event  in  Siberia 
in  1908  appears  to  be  such  an  atmospheric  impact.  Estimates 
of  energy  associated  with  this  event  range  from  3  to  24  megatons 
of  TNT  (51,  52),  powerful  enough  to  produce  an  air  shock  that 
leveled  approximately  80  million  trees  across  2,000  km2  of  for¬ 
est.  At  a  distance  of  60  km,  the  air  shock  was  still  able  to  knock 
down  a  Siberian  trader  (53),  and  thermal  radiation  was  intense 
enough  to  char  his  clothing  (49).  Even  though  the  Tunguska  at¬ 
mospheric  impact  formed  no  known  crater,  it  produced  MSp  (54) 
and  lonsdaleite  (55).  Studies  of  such  atmospheric  impacts  indi¬ 
cate  that  Tunguska-sized  events  up  to  24  megatons  occur  about 
once  every  220  y  (52).  Similar  but  smaller  effects  occurred  during 
the  Trinity  atomic  bomb  test  in  1945,  an  aerial  burst  that  also  left 
no  crater  yet  produced  glassy  surficial  sheet  melt,  along  with 
rounded  and  teardrop-like  glassy  spherules  (56).  Such  an  atmo¬ 
spheric  impact  scenario  is  also  the  best  explanation  for  other  well- 
known  events  with  no  known  craters,  including  the  Libyan  Desert 
glass  field  and  Dakhleh  Oasis  glass  in  Egypt.  In  the  Australasian 
tektite  field  (780  ka),  microspherules  and  tektites  are  strewn 
across  10-30%  of  Earth’s  surface,  producing  the  world’s  largest 
ejecta  field  and  yet,  there  is  no  known  crater.  Wasson  (57)  pro¬ 
posed  that  the  Australasian  field  resulted  from  an  atmospheric 
impact  by  a  comet  approximately  1  km  in  diameter,  striking 
Earth’s  atmosphere  at  an  oblique  angle. 

The  amount  of  kinetic  energy  transferred  during  an  atmo¬ 
spheric  impact  via  air  shock  depends  upon  the  cross-sectional 
area  of  the  object,  its  velocity,  and  its  mass.  Air  shock  pressure 
depends  upon  the  velocity  of  the  object  and  the  density  of  air  at 
altitude.  Shock  front  temperature  is  limited  to  approximately 
20,000  K  by  dissociation  of  air  molecules  ahead  of  front  (58), 
and  effective  duration  of  the  intense  thermal  pulse  can  be  of 
the  order  of  seconds.  Whether  an  object  disintegrates  in  flight 
depends  upon  its  strength,  size,  shape,  velocity,  and  angle  of 
entry.  In  the  case  of  a  comet  that  is  a  dusty  porous  snowball  hav¬ 
ing  little  strength,  a  20-km-diameter  comet  traveling  at  40  km/s 


would  not  disintegrate  in  the  Earth’s  atmosphere;  the  front  of  the 
comet  would  impact  Earth  before  the  shock  from  atmospheric 
impact  reached  the  rear  of  the  comet.  However,  comets  with  di¬ 
mensions  of  tens  of  meters  will  disintegrate  at  high  altitude. 
Weissman  (59)  estimated  that  a  comet  would  have  to  be  >350  m 
in  diameter  to  penetrate  Earth’s  atmosphere  and  form  a  crater, 
depending  upon  angle,  velocity,  etc.  Such  an  event  would  be  at 
least  500  x  more  energetic  than  the  Tunguska  event. 

YD  Impact  Model.  Based  on  current  data,  we  propose  the  following 
preliminary  model  for  formation  of  the  YDB  NDs  and  MSp.  A 
comet  or  asteroid,  possibly  a  previously  fragmented  object  that  was 
once  greater  than  several  hundred  meters  in  diameter,  entered  the 
atmosphere  at  a  relatively  shallow  angle  (>5°  and  <30°).  Thermal 
radiation  from  the  air  shock  reaching  Earth’s  surface  was  intense 
enough  to  pyrolyze  biomass  and  melt  silicate  minerals  below  the 
flight  path  of  the  impactor  (60).  Pyrolytic  products  were  oxidized, 
locally  depleting  the  atmosphere  of  oxygen,  and  within  microse¬ 
conds,  residual  free  carbon  condensed  into  diamond-like  crystal 
structures,  CSp,  carbon  onions,  and  aciniform  soot.  This  involved 
a  CVD-like  process  similar  to  diamond-formation  during  TNT  de¬ 
tonation.  In  some  cases,  carbon  onions  grew  around  the  NDs  and 
other  nanomaterials.  At  the  same  time,  iron-rich  and  silicate  ma¬ 
terials  may  have  melted  to  form  MSp. 

Several  seconds  later,  depending  on  the  height  of  the  thermal 
radiation  source,  the  air  shock  arrived.  NDs,  MSp,  CSp,  and  other 
markers  were  lofted  by  the  shock-heated  air  into  the  upper  atmo¬ 
sphere,  where  prevailing  winds  distributed  them  across  the 
Northern  and  Southern  Hemispheres.  We  suggest  that  the  above 
model  can  account  for  the  observed  YDB  markers. 

Methods 

Core  samples  were  divided  into  multiple  aliquots  for  a  wide  range  of  ana¬ 
lyses.  Details  on  methodology  followed  the  protocol  published  in  Firestone 
et  al.  (2007),  along  with  improvements  as  discussed  in  detail  in  the  manu¬ 
script  and  in  SI  Appendix,  Methods.  These  include  MSp,  magnetic  grains, 
framboidal  spherules,  magnetic  glass,  aciniform  soot,  pollen,  diatoms,  char¬ 
coal,  and  CSps.  NDs  were  extracted  using  the  procedure  published  in  Kennett 
et  al.  (2),  as  further  discussed  in  the  manuscript  and  in  SI  Appendix,  Methods. 
Standard  procedures  were  followed  for  analyses  of  all  proxies. 

Summary 

Synchronous  peaks  in  multiple  YDB  markers  dating  to  12.9  ka 
were  previously  found  at  numerous  sites  across  North  and  South 
America  and  in  Western  Europe.  At  Lake  Cuitzeo,  magnetic  im¬ 
pact  spherules,  CSps,  and  NDs  form  abundance  peaks  within  a 
10  cm  layer  of  sediment  that  dates  to  the  early  part  of  the  YD, 
beginning  at  12.9  ka.  These  peaks  coincide  with  anomalous  envir¬ 
onmental,  geochemical,  and  biotic  changes  evident  at  Lake  Cuit¬ 
zeo  and  in  other  regional  records,  consistent  with  the  occurrence  of 
an  unusual  event.  Analyses  of  YDB  acid-resistant  extracts  using 
STEM,  EDS,  HRTEM,  SAD,  FFT,  EELS,  and  EFTEM  indicate 
that  Lake  Cuitzeo  nanoparticles  are  dominantly  crystalline  carbon 
and  display  d-spacings  that  match  various  ND  allotropes,  including 
lonsdaleite.  These  results  are  consistent  with  reports  of  abundant 
NDs  in  the  YDB  in  North  America  and  Western  Europe. 

Although  the  origin  of  these  YDB  markers  remains  specula¬ 
tive,  any  viable  hypothesis  must  account  for  coeval  abundance 
peaks  in  NDs,  magnetic  impact  spherules,  CSps,  and  charcoal 
in  Lake  Cuitzeo,  along  with  apparently  synchronous  peaks  at 
other  sites,  spanning  a  wide  area  of  Earth’s  surface.  Multiple 
hypotheses  have  been  proposed  to  explain  these  YDB  peaks 
in  markers,  and  all  but  one  can  be  rejected.  For  example,  the  mag¬ 
netic  impact  spherules  and  NDs  cannot  result  from  the  influx  of 
cosmic  material  or  from  any  known  regular  terrestrial  mechan¬ 
ism,  including  wildfires,  volcanism,  anthropogenesis,  or  alterna¬ 
tively,  misidentification  of  proxies.  Currently,  only  one  known 
event,  a  cosmic  impact,  can  explain  the  diverse,  widely  distributed 
assemblage  of  proxies.  In  the  entire  geologic  record,  there  are 
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only  two  known  continent-wide  layers  with  abundance  peaks  in 
NDs,  impact  spherules,  CSps,  and  aciniform  soot,  and  those 
are  the  KPg  impact  boundary  at  65  Ma  and  the  YDB  boundary 
at  12.9  ka. 
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